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Abstract 

Bi-presymplectic chains of one-forms of arbitrary co-rank are considered. The conditions in which 
such chains represent some Liouville integrable systems and the conditions in which there exist related 
bi-Hamiltonian chains of vector fields are presented. In order to derived the construction of bi- 
presymplectic chains, the notions of dual Poisson-presymplectic pair, d-compatibility of presymplectic 
forms and d-compatibility of Poisson bivectors is used. The completely algorithmic construction of 
separation coordinates is demonstrated. It is also proved that Stackel separable systems have bi- 
inverse-Hamiltonian representation, i.e. are represented by bi-presymplectic chains of closed one- 
forms. The co-rank of related structures depends on the explicit form of separation relations. 
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1 Introduction 



The theory of finite dimensional conservative integrable systems has a long history, starting from the 
works of Lagrange, Hamilton and Jacobi in the first half of XIX century. In fact the Hamilton- Jacobi (HJ) 
theory is one of the most powerful methods of integration by quadratures a wide class of systems described 
by nonlinear ordinary differential equations, with a long history as a part of analytical mechanics. The 
theory in question is closely related to the Liouville integrable Hamiltonian systems. The main difficulty 
of the HJ approach is that it demands a distinguished coordinates, so called separation coordinates, in 
order to work effectively. 

There are two efficient and systematic methods of construction of separation variables for dynamical 
systems. The first one bases on Lax representation and r-matrix theory for derivation of separation coor- 
dinates jl] . In this approach the integrals of motion in involution appear as coefficients of characteristic 
equation (spectral curve) of the Lax matrix. This method was successfully applied for separating vari- 
ables in many integrable systems [l]-[6]. The other one is a geometric separability theory on bi-Poissonian 
manifolds [7] -[II], related to the so-called Gel'fand-Zakharevich (GZ) bi-Hamiltonian systems (T51 ITS], 
In this approach the constants of motion are closely related to the so-called separation curve which is 
intimately related to the Stackel separation relations. 

The bi-Poissonian formulation of finite dimensional integrable Hamiltonian systems has been system- 
atically developed for the last two decades (see [UJ and the literature quoted there). It has been found 
that most of the known Liouville integrable finite dimensional systems have more then one Hamiltonian 
representation. Moreover, in the majority of known cases, both Poisson structures of a given flow are 
degenerated. For such systems, related bi-Poissonian (bi-Hamiltonian) commuting vector fields belong to 
one or more bi-Hamiltonian chains starting and terminating with Casimirs of respective Poisson struc- 
tures. The most important aspect of such a construction is its relation to the geometric separability 
theory. Having a bi-Hamiltonian representation of a given system, the sufficient condition for the ex- 
istence of separation coordinates is the reducibility of one of the Poisson structures onto a symplectic 
leaf of the other one. Unfortunately, this procedure is non-algoritmic and has to be considered inde- 
pendently from case to case. Moreover, we do not have a proof that it is always possible for any GZ 
system. Anyway, once the reduction is done, the remaining procedure of the construction of separation 
coordinates is almost algoritmic. The relevance of bi-Hamiltonian formalism in separability theory was 
recently confirmed in |18) . where it was proved that arbitrary Stackel system, defined by an appropriate 
separation relations, has a bi-Hamiltonian extension. 

On the other hand, it is well known from the classical mechanics, that if the Poisson structure is 
nondegenerate, i.e. if the rank of the Poisson tensor is equal to the dimension of a phase space, then 
the phase space becomes a symplectic manifold with a symplectic structure being just the inverse of the 
Poisson structure. In such a case there exists an alternative (dual) description of Hamiltonian vector 
fields in the language of symplectic geometry. So, a natural question arises, whether one can construct 
such a dual picture in the degenerated case, when there is no natural inverse of the Poisson tensor |19j . 
For such tensors the notion of dual prcsymplcctic structures was developed in [20] [21] . 

The prcsymplcctic picture is especially interesting in the case of Liouville integrable systems. As 
was mentioned above, there is well developed bi-Hamiltonian theory of such systems, based on Poisson 
pencils of the GZ type, with polynomial in pencil parameter Casimir functions and related separability 
theory. The important question is whether it is possible to formulate an independent, alternative bi- 
presymplectic (bi-inversc-Hamiltonian in particular) theory of Liouville integrable systems with related 
separability theory and how both theories are related to each other. 

The advantage of formalism presented is as follows. In the bi-Hamiltonian approach the existence of 
bi-Hamiltonian representation of a given flow is a necessary condition of separability but not a sufficient 
one. Contrary, the existence of bi-presymplectic representation of a flow considered is a sufficient condition 
of separability. Moreover, the construction of separation coordinates is a fully algorithmic procedure (in a 
generic case obviously), as the restriction of both presymplectic structures to any leaf of a given foliation 
always exists and is a simple task. For this reason the new formalism presented in the paper seems to be 
relevant for the modern separability theory. 

The present paper develops the general bi-presymplectic theory of Liouville integrable systems when 
the co-rank of presymplectic forms is arbitrary. The whole formalism is based on the notion of d- 
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compatibility of presymplectic forms and d- compatibility of Poisson bivectors. Some elements of that 
formalism was presented in papers |21[ 122] . Here we present a complete picture. Finally it is shown 
that any Stackel system, defined by an appropriate separation relations, has a bi-inverse-Hamiltonian 
representation, what confirms the relevance of presented formalism. 

The paper is organized as follows. In section 2 we give some basic information on Poisson ten- 
sors, presymplectic two-forms, Hamiltonian and inverse Hamiltonian vector fields and dual Poisson- 
presymplectic pairs. In sections 3 the concept of d-compatibility of Poisson bivectors and d-compatibility 
of closed two-forms is developed. Then, in section 4, the main properties of bi-presymplectic chains of 
arbitrary co-rank are investigated. The conditions in which the bi-presymplectic chain is related to some 
Liouville intcgrablc system and the conditions in which the chain is bi-inversc-Hamiltonian are presented. 
Moreover, the conditions in which Hamiltonian vector fields, constructed from a given bi-presymplectic 
chain, constitute a related bi-Hamiltonian chain are also found. In section 5 we prove that arbitrary 
Stackel system, defined by an appropriate set of separation relations, has a bi-inverse-Hamiltonian for- 
mulation. Finally, in section 6, we illustrate presented theory by few representative examples. 

Our treatment in this work is local. Thus, we always restrict our considerations to the domain O of 
manifold M where appropriate functions, vector fields and one-forms never vanish and respective Poisson 
tensors and presymplectic forms are of constant co-rank. In some examples we perform calculations in 
particular local chart from O. 

2 Preliminaries 

Given a manifold M. of dim A'! = m, a Poisson operator H of co-rank r on M. is a bivector n 6 A 2 (A4) 
with vanishing Schouten bracket: 

[n,n] s = o, (2.1) 

whose kernel is spanned by exact one-forms 

kern = Sp{dci}i=i t ... tr . 
In a local coordinate system (x 1 , . . . , x m ) on Ai we have 

while the Poisson property (|2.1[) takes the form 

V (n tf ftir* + ir^n^' + u kl diW l ) = o, ■.= 

— J ox 1 

i 

Let C(M) denote the space of all smooth real-valued functions on M.. A function c £ C{M) is called 
the Casimir function of the Poisson operator n if Hdc = 0. Having a Poisson tensor we can define a 
Hamiltonian vector fields on M. A vector field Xp related to a function F G C(M) by the relation 

X F = UdF, (2.2) 

is called the Hamiltonian vector field with respect to the Poisson operator n. 

A linear combination = Hi + AHo (A S M) of two Poisson operators Ho and Hi is called a Poisson 
pencil if the operator is Poisson for any value of the parameter A. In this case we say that Ho and 
Hi are compatible. When all Casimir functions of are polynomials in parameter A then we say that 
the pencil is of Gel'fand-Zakharevich (GZ) type. 

Further, a presymplectic operator Q on M. is defined by a two-form that is closed, i.e. dfl = 0, 
degenerated in general. In the local coordinate system (x 1 , . . . , x m ) on M. we can represent ft as 

rn 

fl = flijdx 1 A dx J , 

i<j 



2 



where the closeness condition takes the form 

dittjk + dk^ij + djtt k i = 0. 

Moreover, the kernel of any presymplectic form is an intcgrable distribution. A vector field X F related 
to a function F £ C(M) by the relation 

ilX F = dF (2.3) 

is called the inverse Hamiltonian vector field with respect to the presymplectic operator O. 

As in the case of presymplectic forms their linear combination is always presymplectic, hence the 
notion of compatibility, as it was defined for Poisson tensors, does not make sense. We will come back to 
this problem in the next section. 

Any non-degenerate closed two form on M. is called a symplectic form. The inverse of a symplcctic 
form is an implectic operator, i.e. invcrtiblc Poisson tensor on M. and vice versa. 

Definition 1 A pair (n, 0) is called dual implectic- symplectic pair on M. if II is non- degenerate Poisson 
tensor, fl is non-degenerate closed two-form and the following partition of unity holds on TAi, respectively 
on T*M: I = 110. and I = On. 

So, in the non-degenerate case, dual implectic-symplectic pair is a pair of mutually inverse operators 
on M.. Moreover, the Hamiltonian and the inverse Hamiltonian representations are equivalent as for any 
implectic bivector n there is a unique dual symplectic form SI = n _1 and hence a vector field Hamiltonian 
with respect to n is an inverse Hamiltonian with respect to Q. 

Let us extend these considerations onto a degenerate case. In order to do it let us introduce the 
concept of dual pair as it was done in |20j . Consider a manifold M. of an arbitrary dimension m. 

Definition 2 [20] A pair of tensor fields (n, fi) on M. of co-rank r, where II is a Poisson tensor and f2 
is a closed two-form, is called a dual pair (Poisson-presymplectic pair) if there exists r one-forms dci and 
r linearly independent vector fields Zi, such that the following conditions are satisfied: 

1. kern = Sp{dci : i = 1, ...r}. 

2. kerJl = Sp{Z % : i = l,...r}. 

3. Zi(cj) = Sij , i = 1,2 . . . r. 

4- The following partition of unity holds on TA4, respectively on T*A4 

r r 

i = un + ^2z l ®dc l , i = on + ^dcj® z u (2.4) 

i=l i=l 

where <S> denotes the tensor product. 

A presymplectic form plays the role of an 'inverse' of Poisson bivector n in the sense that on any 
symplectic leaf of the foliation defined by kern, the restrictions of il and n are inverses of each other. 
More information on geometric interpretation of dual pairs the reader can find in [20j . Contrary to the 
non-degenerated case, for a given Poisson tensor n the choice of its dual is not unique. Also for a given 
presymplectic form f2 the choice of dual Poisson tensor is not unique. We will come back to that problem 
at the end of this section. 

For the degenerate case the Hamiltonian and the inverse Hamiltonian vector fields are defined in the 
same way as for the non-degenerate case, but for degenerate structures the notion of Hamiltonian and 
inverse Hamiltonian vector fields do not coincide. For any degenerate dual pair it is possible to find a 
Hamiltonian vector field that is not inverse Hamiltonian and an inverse Hamiltonian vector field that 
is not Hamiltonian. Actually, assume that (n, fi) is a dual pair, Xp = HdF is a Hamiltonian vector 
field and dF = CIX F is an inverse Hamiltonian one-form, where X F is an inverse Hamiltonian vector 
field. Having applied ft to both sides of Hamiltonian vector field, n to both sides of inverse Hamiltonian 
one-form and using the decomposition (|5.3[) we get 

r r 

dF = n(X F )+Y J Z l (F)dc t7 X F ^ X F -Y J X F {c,)Z l . (2.5) 

i=l i=l 
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It means that an inverse Hamiltonian vector field X F is simultaneously a Hamiltonian vector field Xp, 
i.e. X F = Xp, if dF is annihilated by ker(fi) and X F is annihilated by ker(II). Moreover, for any dual 
pair (II, fi), the following important relations hold [20] 

[z it z s ] = o, l Xf u = o, L Zi u = o, L XF n = o, L Zz n = o, (2.6) 

for i,j = 1, r, where is the Lie-derivative operator in the direction of vector field X and [. , .] is a 
commutator. 

Let us return to a 'gauge freedom' for a duality property. In other words: given a dual pair (II, fi) 
how can we deform fi to a new presymplectic form fi' so that (II, fi') is again a dual pair, or how can we 
deform II to a new Poisson operator II' so that (II', O) is also a dual pair? 

Lemma 3 i2Uf Let II be a fixed Poisson tensor and ft be a dual presympectic form. Assume that dci £ 
kcrll, Zi € kerfi and Zi{c.j) = 5{j. Define 

fi' =Q, + Y 1 df i Adc i , 

i 

where fi £ C(M). Then (II, Q) is a dual pair, with ker(fi') = Sp{Z i = Zi — Udfi}, provided that 

Z i (f j )-Z j (f i )+U(df i ,df j ) = for all (2.7) 

Lemma 4 \20f Let O be a fixed presymplectic form and II be a dual Poisson tensor. Assume that 
Zi 6 kerfL dci € kerll and Zi(cj) = dij. Define 

n' = n + ^A^, (2.8) 

i 

where Ki are vector fields such that 

Ki=ILdFi, dF t = nK t Z j (F i ) = 0, K J (c l )=0, i,j = l,...,r, (2.9) 

for some functions Fi £ C{M). Then, (LT',0) is a dual pair, with ker(LV) = Sp{dc[}, c[ = Cj + Fi, 
provided that 

Ct(K i ,K j )=0 forallij. 

Poisson tensor II, considered as the mapping II : T* M. — > TA4, induces a Lie bracket on the space 
C{M) 

{.,.} n :C(M) xC(M) ^C{M), {F,G} n = (dF,ndG)=n(dF,dG), (2.10) 

(where (., .) is the dual map between TM. and T*A4) which is skew-symmetric and satisfies Jacobi 
identity. It is called a Poisson bracket. Jacobi identity for ()2.10[) follows from the property (|2.Ij) of II. 

When a Poisson operator II is nondegenerate its dual ft is its inverse ft = II -1 . Moreover, any 
Hamiltonian vector field with respect to II is simultaneously the inverse Hamiltonian with respect to Q 
and Xp = X F . Hence, a symplectic operator ft defines the same Poisson bracket as the related Poisson 
operator n 

{F, G} n d = n(X F , X G ) =< ftX F , X G >=< dF, X G >=< dF, UdG > (2. II) 

= U(dF,dG) = {F,G} n . 

What is important, when II is a degenerate Poisson tensor and f2 is its an arbitrary dual two-form, 
the formula (|2.11j) is still valid. It follows from the fact that although Xp ^ X F , but < VtXp, UdG >=< 

nx F ,ndG >. 

Finally, we remind the reader two identities important for further considerations. Let n be a Poisson 
bivector and be a closed two-form, then 

£n 7 n + nd 7 n = o, L T o. = d{Q,T), (2.12) 

where r G TM and 7 £ T*M. 
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3 D-compatibility of closed two-forms and Poisson bivectors 

In the following section we develope a concept of d-compatibility which is crutial for our further consid- 
erations. Let us start with a non degenerate case. 

Definition 5 We say that a closed two-form fli is d- compatible with a symplectic form f^o if Ho^IiHq is 
a Poisson tensor and Ilo = fip 1 is dual to Qq. 

Definition 6 We say that a Poisson tensor Hi is d-compatible with an implectic tensor n if £IqTIi£Iq 
is closed and fl = IIq~ 1 is dual to TLq. 

Now, the following lemma relates d-compatible Poisson structures, of which one is implectic, and 
d-compatible closed two-forms, of which one is symplectic. 

Lemma 7 fMI 

Let (Ilojfio) be a dual implectic- symplectic pair. 

(i) Let a Poisson tensorlli be d-compatible with Hq. Then fio and ^i = ^ollifio ore d-compatible closed 
two-forms. 

(ii) Let a closed two-form fli be d-compatible with S7o . Then Hq and Hi = Ilofiillo are d-compatible 
Poisson tensors. 

Let us extend the notion of d-compatibility onto the degenerate case. 

Definition 8 A closed two-form tti is d-compatible with a closed two-form fio if there exists a Poisson 
tensor Hq, dual to fioj such that Tloftillo is Poisson. Then we say that the pair (f^o , Oi ) is d-compatible 
with respect to Hq. 

Definition 9 A Poisson tensorlli is d-compatible with a Poisson tensor Hq if there exists a presymplectic 
form fio, dual to Hq, such that ^ollifio is closed. Then we say that the pair (Ilo, III) is d-compatible 
with respect to £Iq. 

Comparing the notions of compatibility and d-compatibility for Poisson pair (IIo,IIi) we will show 
that when Ilo is non-degenerated both notions are equivalent, but for a degenerate case the notion of 
d-compatibility is the stronger one. Actually, let us consider the following identity, proved in [21| . 

£(n 1+ Anoh( n i + An o) + (Hi + An )d 7 (ni + An ) 

= A{L T (ftoIIiOo) - d(O nifi r) - EJ^o(£z i n 1 )f2 ]r A dc t (3.1) 

-E^cOno^nOfio}, 

where IIo,IIi are Poisson tensors, (IIo,rio) is a dual pair, where dci £ kerLTo, Zi £ kerfioj T £ TA4 and 
7 = SloT £ T*A4. Assume first that Ilo an d LTi are d-compatible with respect to ilo- Then iloLTifio is 
closed and 

L T (Q nifio) -d(f2 niQ T) = 0, teTM. (3.2) 
In particular, for r = Zi, relation (|3.2j) gives 

n (L Zi Hi)n = o, i = i,...,r. (3.3) 

Hence 

£(n 1+ An b( n i + W + (Hi + An )d 7 (ni + An ) = (3.4) 

and LTi + Alio is Poisson. On the other hand, from the compatibility relation (|3.4I) the d-compatibility 
(|3.2p follows under additional conditions (|3.3j) . 



5 



(k) 

Theorem 10 Let a Poisson tensor Hq and a closed two-form CIq form a dual pair, where Y G ker 17q 



dH^ 1 e kcrn and Y {k) (H^ m) ) = S km , k,m=l,...,r. 

(i) If Hi is a Poisson tensor d-compatible with Hq with respect to 17o, then forms 17o an d ^1 = 17onil7o 
are d-compatible with respect to Hq. 

(ii) If fii is a closed two-form d-compatible with 17o with respect to Hq, then Poisson tensors Hq and 
Hi = nofiillo are d-compatible with respect to 17q, provided that 



H niY Q {k) = H dF {k \ k = l,...,r 
for some functions F^ k > £ C(M) and 

17 1 (Y (fe) ,y o (m) ) + Y {k) (F^) - Y {m) (F^) = const, 



k, m = 1, 



(3.5) 



(3.6) 



Proof. 

(i) fii is closed as III is d-compatible with IIo. Then, nol7i±Io = nol7onil7oIIo is Poisson (as was shown 
in [H]). 



(ii) From the d-compatibility of S7q and 17 1 it follows that III is Poisson. Then, 



17QIT170 = 17 n 17inof2o = (J - E fe ^ fe) ® Y^Q^I - E,„ C ® dH t 



,(fe)N 



r(m) 



r( m A 



iH^AniiY^ + l^Sli 



r (m) v (fe)\ 



r(fc) 



(m) 



From the assumption (|3.5[) and decompositions ()5.3|) it follows that 



i7iF w = ^ fcJ +E m i7i(Y (fc) ,y (m) )-^ m; (^ w ) 



(m) 



hence, 



I7 riii7 = S7! + J2k dH o k) A dF 



Efc,; 



= i7 x + Efcdffo A dF(fe) 
+ jEm. [ox(y (fe) ,F (m) ) - y (m) (F( fc )) + yftfH) 



and under condition (|3 .€>[) 17oriil7o is closed. □ 

The important for further considerations special case occurs when 



i7i (yf»F H ) = o, r«(F (ro) ) = F (m) (^ (fe) )> = l,..,r. 



(3.7) 



Theorem 11 Le£ a Poisson tensor Hq and a closed two-form 17q form a dual pair, where Y £ kcr 17q 



dH^ ekern and Y (K > (H^" 1 ') 



k, m = 1. 



(i)If Hi is a Poisson tensor d-compatible with Hq with respect to 17q and 



X {k) = HidH (k) = H dH\ 



r(fc) 



fc = 1, 



(3.8) 



are bi-Hamiltonian vector fields for some functions H[ k \ then 17o and 17i = 17oIIil7o + Efe'^i ^ dH^ 
is d-compatible pair of presymplectic forms with respect to Hq. 



(ii) //17i is a presymplectic form d-compatible with 17o with respect to Hq and 



(3.9) 



are bi-presymplectic one-forms, then Poisson tensors Hq and Hi = noOiIIo + 5Zfc-^^ A i j are ^~ 
compatible with respect to ilo if 

X [k) = n Oiy (fe) = n dF (fe) , noOiF^ 5 = U Q dG (k) , k = 1, ...,r, (3.10) 
/or some functions F^ k \G^ £ C(jVl) and 

no( yW )F M ) = 0) fc, m =l,..., r , (3.11) 

F (m) (F (fc) ) = F (fe) (F (m) ), fc,m = l,..,r, (3.12) 

y 1 (m) (fff ) )=y {i,) (FW) 1 fc,m = l,...,r. (3.13) 

Proof. 

(i) fl\ is closed as III is d-compatible with n . Then, n f2in = no^oriifiono is Poisson (as was shown 
in[H]). 

(ii) From f|3.10[) we have 

n 1 (Y o (k) ,Y o {m) ) = (3.14) 

and by previous theorem part (ii) the form f2oIIif2o = fiono^inoflo is closed under condition ()3.12[> . 
Moreover, (|3.14p yields 

n Q Y^ = dF« _ J2 m Y^(F^)dH^ m \ 

fiilf = dG^ - T m Yt\G (k) )dH^ n) . 
Conditions p. lip and (|3.13p arc sufficient for III to be a Poisson tensor. From p. lip it follows that 

Qo[Y^ k) ,Y^ m) ) =0=> HoidF^ldF^) = =*> [X {k \X^ m) ] = 0. (3.15) 

Now we show that the Schoutcn bracket of Hi is zero. As n f2iIIo is Poisson (it follows from compatibility 
of Qq and fix), we have 

pi, ni] s = 2E fc [n f7 1 n ,x( fe ) a y^% + E*,J* (fc) a F„ (t, ,i (m) a y (m) ] s , 
[n Oin ,xW a y { % = Y^ k) AUodi^x^no-x^ AUodin^^Uo, 

[X<V Ay W ,I W AY im) ] s = 2X^AY (m) A [Y (k \x^]. 
In last equality we used the fact that [l^*^, ] = and relation (|3.15p . Now, 

= [yW^ofii^] = L y(fc) (n ^ 1 )y (m) = n (L y(fc) o 1 )y (m) 
= n d(f7 lM y (fe) )y (m) = n (d/^>)y (m) 



Zid(YW(FV))AdH^ 



r 



= n 

= -n d(y (fc) (^ M ))- 

From 

yW = X (fe) + Y, m Y^ k) {H^ m) )Y^ m) 

we have 
and 

fiiXW = dG^ - J2 t Y Q {k) (G^)dH { a l) - J2 t Y} k) (H^)dF^ + Zi,jY{ k) (H®)Y® (F^)dH^ . 
Hence, 

n d(fiiX«)n = -no [Em<*(*i fc) (#o m) )) a dF^~ 
= E m no« ( Vr ) ))A^ w 

and then [ni,IIi]s = under condition (|3.13[) ,D 



n (: 
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4 Bi-presymplectic chains 

Now we are ready to investigate main properties of bi-presymplectic chains. 

Theorem 12 Assume we have a pair of presymplectic forms (f2o, Oi), d-compatible with respect to some 
Ho dual to Qo, both of rank 2n and co-rank r on O C A4. Assume further, that they form bi-presymplectic 
chains of one-forms 

p\ k > = n Y} k) = n.Y^l i=l,2,...,n k (4.1) 

where k = 1, ...,r, n\ + ... + n r = n and each chain starts with a kernel vector field Y" of O and 
terminates with a kernel vector field Yn k of f2i . Then 

(i) 

n (Yf k) ,Y} m) ) = fi!^,^') = 0, (4.2) 

for k,m= 1, ...,r, i = 1,2,..., nk, j = 1,2,..., n m . 
Moreover, let us assume that 

x\ k) = n (3 {k) = n dH$ k \ (4.3) 
for k = 1, r, i — 1, 2, . . . , rik, which implies 

fjto = dffW - ^(fll*))^, (4.4) 

yW = jr« + E m i r * w (^ m) )i'o <m) . (4-5) 

where HodHo = 0. Then, 
(ii) 

U {dH[ h) , dH f° ) = 0, [Xf 5 , Xj m) ] = (4.6) 
and equations (RTTl) define a Liouville integrable system. 



Additionally, if 
and 



Y^\H ( T ) )=Y^ n) {H[ k) ) (4.7) 



^(ffW) = yW(^)), (4.8) 
where k, m = 1, r, i = 1, 2, . . . , n m , then 

(Hi) vector fields x!- form bi-Hamiltonian chains 

X\ k) = U dHl k) = n x dH\%, i = 1,2, . . . , n (4.9) 

where 

n x = noOxiio + E m x[ m) a y (m) , (4.io) 

k,m = 1, r, i = 1,2, . . . ,ftfc and ni + ... + n r = n. The chain starts with H^ k \ a Casimir of Hq, and 
terminates with , a Casimir of Hi. Moreover the Poisson pair (Iio,IIi) is d-compatible with respect 
to VLq. 

Proof. 

(i) From (|4.1[) we have 

n o(F w, F .(-)) = n ( F w, F w )) 

1 (F.( fe ),F.(-)) = n 1 (F«,y.M) 
f2o( yW ) vM ) = O l( yW ) y(-) ) . 
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Then, (g21) follows from 

f2o (yW ! y/-)) = o ! fli(yW^ n) ) = o. 

(ii) From properties of dual pair (IIo, f2o), if X\ = TlgdH^ then 

n (dff w , dH f°) = n (x^ , ) . 

On the other hand as x\ k) = Y % (k) + E m ot$i}Yo m \ where a^} are an appropriate functions, it follows 
that 



no(xl k \xl m) ) = n (Y z {k \Yf m) ). 



(iii) We have 



xj fe) = HndH$ k) 



of 

= nofiiy S = noouxw + E ro ^ ( -l(^ m) )io (m) ) 
= nofiiiiodfl^* 5 ! + E m ^ ( _i(^ m) )^ m) ) 

(n o ino + E m ^ (m) a r Cra) 
= nidfl 1 ( * ) 1 . 

Moreover, IIo and III are d-compatible Poisson tensors provided that (|4.7|) is fulfilled. We also have 

<P n o x (yw - Eji?^ 1 )^) + E m ^, (ro) (^ ) )^ m) 

□ 

Notice, that in a special case, when 

F (fe) (ffH) =0, (4.11) 

for all admissible values of fc, m and i, chains (|4.ip are bi-inverse-Hamiltonian as = dH^ k \ Obviouslu 
X\ k) are not bi-Hamiltonian until Y^ k) {H^ n) ) ^ 0. 

Finally we show that arbitrary Liouville integrable system which has a bi-prcsymplcctic representation 
on (2n + r)-dimcnsional phase space, has also quasi-bi-Hamiltonian representation on any symplectic leaf 
of its Hamiltonian structure IIo ■ Actually, from (|4.1I) , (14. 3p and (|4.4[) follows that 

U dH^ = ILA (y^\ + E m yt\H^)dHt ] ) 

(k) ^ f v (k) , rrH\ n v (m) v (m)frT(k)^,rr(m) 



= no [n.x^i + E m (^i«"0^i^ mj + >TWW 
= noiiiUodHjp! + E m Y ^l( H o m) WH[ m \ 

hence on (2n + r)-dimcnsional phase space we have quasi-bi-Hamiltonian representation 

r 

li^HfX = IL dH$ k) + J2 F^Il dHi m) , (4.12) 

m— 1 

where 

n a = no^rno, F^ m) = -y/ fe) (^ m) ). 

Notice that both Poisson structures IIo and Hi = Hofiillo share the same Casimirs Hq , so the 
quasi-bi-Hamiltonian dynamics can be restricted immediately to any common leaf of dimension 2n 



■Kxdhf\ = ir dh {k) + J2 F^ n) TT dh[ m \ i = 1, n, (4.13) 

m— 1 
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where iTi and h\ are restrictions of 11; and H\ , respectively. Hence we deal with a Stackel system whose 
separation coordinates are eigenvalues of the recursion operator N = t^i^q 1 |23] . provided that N has n 
distinct and functionally independent eigenvalues at any point of O C Ai, i.e. we are in a generic case. 
We will come back to separable systems in next sections. 

The advantage of bi-presymplectic representation of Liouville integrable system, when compared to 
bi-Hamiltonian ones, is that the existence of the first guarantees that the system is separable and the 
construction of separation coordinates is purely algorithmic (in a generic case) , while the bi-Hamiltonian 
representation does not guarantee the existence of quasi-bi-Hamiltonian representation and hence sep- 
arability of the system in question. Moreover, the projection of the second Poisson structure onto the 
symplectic foliation of the first one, in order to construct a quasi-bi-Hamiltonian representation, necessary 
for separability, is far from being trivial non-algorithmic procedure |24| . 



5 Separable Stackel systems 

Consider a Liouville integrable system on a 2n-dimensional phase space M, i.e. a set of hi G C(M), 
i = l,...,n which are in involution with respect to some Poisson tensor 7r , and related n Hamiltonian 
dynamic systems 

u ti — 7r dhi = Xhi, i=l,...,n, (5-1) 

where u G M and Xh t are respective Hamiltonian vector fields. 

The Hamilton- Jacobi (HJ) method for solving (|5.1j) essentially amounts to the linearization of the 
latter via a canonical transformation (when u G O C M. is some local canonical chart) 

u=(q,p)^ (b, a), a,i = hi, i = l,...,n. (5.2) 



In order to find the conjugate coordinates b l it is necessary to construct a generating function W(q, a) of 
the transformation (|5.2[) such that 

The function W(q, a) is a complete integral of the associated Hamilton- Jacobi equations 



h i {q\...,q n ,^,...,^)=a i , i = l,...,n. (5.3) 
In the (b, a) representation the ^-dynamics is trivial: 

(a,j) ti = 0, {b 3 ) u = Sij, 

whence 

V(Q,a) = i^=ti+7i, j = 1, . . . , n, (5.4) 

where jj are arbitrary constants. 

Equations (|5.4j) provide implicit solutions for (|5.1[) . Solving them for qj is known as the inverse Jacobi 
problem. The reconstruction in explicit form of trajectories q J = q J (fa) is itself a highly nontrivial problem 
from algebraic geometry. 

To solve the system (|5.3|) for W in a given local coordinate system is a hopeless task, as (|5.3|) is a 
system of nonlinear coupled partial differential equations. In essence, the only hitherto known way of 
overcoming this difficulty is to find distinguished canonical coordinates, denoted here by (A, /z), for which 
there exist n relations 

ip i (\ l ,[j, i ;ai,...,a n ) = 0, i = 1, . . . ,n, o, € K, det f^- ^ 0, (5.5) 

such that each of these relations involves only a single pair of canonical coordinates pQ . The determinant 
condition in (|5.5[) means that we can solve the equations (|5.5[) for a; and express a; in the form a, = 
hi(\,fi), i = l,...,n. 
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If the functions Wi(X % ,a) are solutions of a system of n decoupled ODEs obtained from (|5.5[) by 
substituting fa = dW '}^i 

<Pi [X l ,fa = dWi jxi ' a \ ai,...,a n ^j = 0, i = l,...,n, (5.6) 

then the function 



W(\,a) = J2. Wi(\\a) 
* — »&= 1 



is an additively separable solution of all the equations (|5.6[) . and simultaneously it is a solution of all 
Hamilton- Jacobi equations (|5.3[) . The distinguished coordinates (A, fx) for which the original Hamilton- 
Jacobi equations (|5.3[) arc equivalent to a set of separation relations (|5.6p arc called the separation 
coordinates. 

In what follows we restrict ourselves to considering a special case of (|5.5[) when all separation relations 
arc affine in ft^: 

n 

J2sHx l ,fa)h k ^^(X\fa), i = l,...,n, (5.7) 

fc=i 

where and ^ are arbitrary smooth functions of their arguments. The relations (|5.7[) are called the 
generalized Stackel separation relations and the related dynamical systems are called the Stackel separable 
ones. The matrix S — (5?) will be called a generalized Stackel matrix. To recover the explicit Stackel 
form of the Hamiltonians it suffices to solve the linear system (|5.7p with respect to hi. 

If in (|5.7p we further have S^{X l ,fa) = S {X l ,fa) and ipi{X l ,fa) = ip(X l ,fa) then the separation 
conditions can be represented by n copies of the curve 

n 

^Ts^A,^ = <A(A, M ) (5.8) 
fc=i 

in (A, n) plane, called a separation curve. The copies in question are obtained by setting A = A 1 and 
/x = fa for i = 1, ... ,n. 

For further convenience, let us collect the terms from the l.h.s. of (|5.7[) as follows: 



Y,d(X\fa)h^(X l )=^ t (X\fa), i = l,...,n, (5.9) 



k=l 

where 



and impose the normalization ^(A z ,/i;) = 1. 

Some informations about the classification of Stackel systems (|5.9[) the reader can find in Q~8] . 



6 Bi-inverse-Hamiltonian representation of Stackel systems 

As recently proved in [T5], the Stackel Hamiltonians defined by separation relations (|5.9[) admit on M 
the following quasi bi-Hamiltonian chains in (A, /x) representation 

r 

7ridh\ k) =-K dh ( Q l +Y J Fi k ' l) ^dh < (\ h ( £ +1 = 0, k = l,...,r, i = l,...,n k , (6.1) 
where itq is a canonical Poisson tensor 

d d 
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iv i is a noncanonical Poisson tensor of the form 

— A — 



(k I) 

compatible with 7ro, and the expantion cefficients ' are solutions of the set of linear algebraic equations 

r 

^^(A^ W )F (M) (A^) = ^.(A^ W )(A j r ! , 3 = !,•••, n, J = l,...,r, (6.2) 
fe=i 

where 

n* 

i=l 

Let us consider the following symplectic forms on M 



F W) (A) = ^ A n *-*if ,0 , m + • • • + n r = n. 



w = - ^ dA* A dfii, wi = — A l dA' A d/ij. 

Observe that (7To,wo) constituts non degenerate dual implectic-symplectic pair as u>q = 7Tq~ , no and 
7Ti = 7ToWi7To are d-compatible with respect to wo and u/q and cji = cjotti^o are d-compatible with 
respect to ttq. Besides, quasi bi-Hamiltonian chains (|6.1[) have equivalent quasi bi- inverse- Hamiltonian 
representations 



W if. wo ^) +1 = 0, fc = l,...,r, i = l,...,n fc) (6.3) 



wix^ = wo a;. 



where 

(*) jlW Jt(fc) (fe) 

x • = 7r an ■ , an ■ = w :c ■ . 

Let us lift the whole construction to the extended phase space M — > J\4 : (A, /x) — > (X,fx,c), where 
dim M. = 2n + r. Then, on A4 : ljq — >• fioj - > IIo, both degenerated, where 

kerf2 = Wcf }, fc = l,...,r, F (,s) = — , n Y (fc) = 

Ccfc 

and 

kerllo = Sp{dc k }, k=l,...,r, Tl dc k = 0, r (fc) (c,-) = ^. 
Obviously, (IIo, ^o) is a dual Poisson-presymplectic pair on JW. In the same fasion we lift 

wi -> fiiu, 7Ti — » IIiu, scf ) Af \ 

where keitlm = kerfio and kerlliD = kerllo. On M. quasi bi-inverse-Hamiltonian chains (|6.3|l take the 
form 

r 

ftiDAf =Qo^ ) 1 + 5Z^ (fe ' °o^i (0 , fc = l,...,r, i = l,...,n k . (6.4) 
Let us define the following prcsymplectic two-form 



fi x = fi^ + ^ A dc fc (6.5) 



fe=i 

and the set of vector fields 

rW=xf)-^FW)y«. (6.6) 
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Then, wc have 



n Yg\ = dh[% 



1=1 

(n t J2 dh i A dc ^ k) + E - E Fi u) dh { ! ] 
1=1 1=1 i=i 



;=i ;=i (=i 

j^F^dh^ + £ F^ m \^\h^)d Cl -J2F^dh^ 

1 = 1 l,m=l 1=1 

r(k) 



a.s 



= fill? 

fii^o = Erf A dc ^ Y o l) = dh i^ 

k=l 

F «(^) )=0) f/%) = -if ^ Yt\c k ) = 5 mk . 
Hence, on A4, differentials dh[ form a bi-invcrsc-Hamiltonian hierarchies 

QoY™ =dh\% =n 1 Y l (k \ i = 0,l,2,...,n fc , fc = l,...,r, (6.7) 

which starts with a kernel vector field Yg*^ of flo and terminates with a kernel vector field Ynf of fii. 
Indeed 

QxYW = (flu, + £ dfc<™> A - £ ^F M ) 

m— 1 m— 1 

= £ - E F&^dh^ = 0. 

m— 1 m— 1 

Moreover, flo and f2i are d-compatible with respect to Hq, as 

iL^no = n o 1Z5 n = n 1D 

which is Poisson. According to theorem [T2l vector fields are not bi-Hamiltonian as Y^ k \h^) = 

_p(k,l) ^ Q 

In order to construct on Ai bi-Hamiltonian representation of considered Stackel systems, one has to 
extend the original Hamiltonians 

h W ^ H (k) = ^^F^d, i = l,...,n. (6.8) 

i=i 

Then, on A4, vector fields 

r 

K (k) = x (k) _ Uod( ^2 F^a) (6.9) 
i=i 

form a bi-Hamiltonian chains 

n dHl%=K < ^ 1 = n 1 dH\ k \ i = 0,l,...,n fe , k = l,...,r, (6.10) 
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whore; 

IIi = n 1D + J2 k { i n) A (6.11) 

m— 1 

is a Poisson tensor compatible with IIo one. Each chain starts with the Casimir of IIo, i.e. Hq = c^, 
and terminates with the Casimir of IIi, i.e. Hn k . The details of the construction the reader finds in [18] . 
Poisson tensors IIo and III ar e d-compatiblc with respect to 57o as 

SloIIifJo = r2olli£)0o = fiizj 

is closed. As was proved in [21] , bi-Hamiltonian chains (|6.10[) have no bi-presymplcctic counterparts as 
the conditions (|4.7[) are not satisfied (see also theorem ITT]) . Indeed 

The only exception is the case of co-rank one (r = 1), as then (|4.7j) is trivially fulfilled. 

7 Examples 

Here we illustrate the presented theory with three examples of separable systems, each of three degrees 
of freedom. Two of them are classical Stackel systems with separation relations being quadratic in mo- 
menta, while the third example has separation relations cubic in momenta. 
Example 1. 

Consider the separation relations on a six-dimensional phase space M given by the following bare 
(potential-free) separation curve 

h x \ 2 + h 2 \ + h 3 = \p 2 ■ 

This curve corresponds to geodesic motion for a classical Stackel system (of Benenti type [21]). As in this 
example k = 1, we use the notation = hi. The transformation (A, /i) — > (q,p) to the flat coordinates 
of associated metric follows from the point transformation 

- qi = -A 1 - A 2 - A 3 , 
02(g) = kql + q 2 =A 1 A 2 + A 1 A 3 +A 2 A 3 , 
03(g) = |gi92 + ?3 = -A 1 A 2 A 3 . 

In the flat coordinates the Hamiltonians take the form 

Ei =PiP3 + \p\, 

E 2 = p\P2 + \q\p\ + \q\pipz - \q2P2p3 - \q-ip\-, 
E 3 = \p\ + yipj + \q\p\ + \q\P\Pi + \qiP\Pz 

- (l<7l<?2 + g 3 )P2P3, 

and admit a quasi bi-inversc-Hamiltonian representation (16. 3[) 

= uj x 2 + F1W0X1, 
UJ1X2 = + F 2 w xi, 

LOlX 3 = F 3 W Xl, 

with the operators loq and uj\ of the form 

-0 = -o- 1 = ( I ~q 3 ) , (7.1) 
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hP2 


1P3 


^91 


\<l2 


</:i 


\ 









"0 


-1 





J* 






-\vz 











-1 








-o<?l 


1 


















-2 12 





1 













I 


-?3 


-592 













/ 



where 4 is an 3 x 3 unit matrix, the expansion coefficients f\ X ^ = Ff. 



(7.2) 



Ft 



-io 2 

4^1 



<?2, 



-F3 = -\qiqi 



?3 



and Hamiltonian vector fields = itodhi, i = 1, 2, 3. 

On the extended phase space M of dimension seven, with an additional coordinate c, the differentials 
dhi form a bi-invcrsc-Hamiltonian chain 

ft F = 

fi Yi = dhi = QiY 
n Y 2 = dh 2 = fiin 

fi y 3 = * 3 = o x y 2 
o = Oir 3 , 

with prcsymplcctic forms 



CJ 














dhi 








d-compatible with respect to ITo and vector fields 

lo = (0,...,0,l) T , Y i = X i -F i Y , i= 1,2,3, 

where = Ilod/ii. 
Example 2. 

Consider the separation relations on a six-dimensional phase space given by the following bare separation 
curve 

X 2 (h^X + h^) + h^ = ^ 2 

representing geodesic motion for a classical Stackcl system (this time of non-Bcnenti type [25]). Using 
the coordinates, the Hamiltonians, and the functions Ui from the previous example we find that 

^ = -±f», 
h { } ] = hi - ^h 2 , 

2 1 CT 2 Z ' 

h{ 2) =h 3 - —h 2 

1 ° a 2 



(k) 

and thus we see that the Hamiltonians h\ ' are related to hj through the so-called generalized Stackcl 
transform (see |26| for further details on the latter). They admit a quasi bi- inverse- Hamiltonian repre- 
sentation (16.31) 



(i) 

UJiX\ 

(i) 

L>JlX 2 

(2) 
LOlX\ 



(1) 
U!qX 2 

#(1,1), 

F 



(2,1), 



F 1 

,(*) 



(1,1) 



UJ X\ 



(1) 



F 



(1,2) 



F 



(1,2), 



1 

(2) 



LO X\ 



(2) 



Z7(2,2) (2) 

F; 'ojqx) 



1 wo^i t x 1 WU^l 
with the presymplectic forms (|7.1[) , (|7.2p . the expansion coefficients 



*i 



(1,1) 



-(7i 



£3 



(1,1) 



-<7 2 



(T2 



(2,1) 
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and Hamiltonian vector fields = irodh\ 

On the extended phase space A4 of dimension eight, with an additional coordinates c\ and c 2 , the 
differentials dh^ form a bi-inverse-Hamiltonian chains 



,(fe) 





= 






= dhP 


= ^„ (1) 




= dh 2 1] 


= o.yw 










o r (2) = o 

= fixF/ 



(2) 



with the presymplectic forms 





( ^0 





O = 













Qi = 



dh[ 1] dh\ 



(2) 







d-compatiblc with respect to Ilo and vector fields 



Y, 



(i) 



(o,...,o,i,o) T , r (2) = (o,...,o,o,i) 



^(M)y(i) _ jT^i^y 



(l,2)v(2) 



*2 — A 2 - ^2 



1 n 



P (l,2) v (2) 



.(2) 



(2,l) v (l) 



,(2,2) v (2) 



Example 3. 

Consider separation relations on a six-dimensional phase space given by the following bare separation 
curve, cubic in momenta, 

The transformation (A,/i) — > (q,p) to new canonical coordinates in which all Hamiltonians are of a 
polynomial form is obtained from the following two transformations: 

Ui = 3q 2 -3g 3 , 

"2 = -qiP2 - qiPs, + 3? 2 + 5qf - 6q 2 q 3 , 
u-3 = ~ql ~ 9qfq3 + qiq3P2 + qiq3Ps - 6qiQ2 
+qlpi +3q 2 q%, 



and 



"2 



1 

3q 2 - 2g 3 
9i 



t* = P3 + f P 2-^ + 3^-4g 2 , 

qi qi 



U\ = A1+A2 + A3, 

u 2 = A1A2 + A1A3 + A2A 3 , 

U3 = A1A2A3, 

Hi = vi\f + v 2 \i + v 3l i = l,2,3. 
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hi 1 



In the (g,p)-coordinates the Hamiltonians take the form 

ill 1 O O O O /1 

_ P2R3 + 3P2 + P3 _ 7 9iP3 - - 3q 2 pi + lSqiq 2 + 13q 1 + I2q 3 qiq 2 

\2q\q 2 + 8qfq 3 - 2q\p x + {-&qiq 2 - 4<?i q 3 )p 3 + pip 3 , 

1 + \v\p-i + §fP2 ~ I1P3 ~ iliPl ~ Q2P1P2 
+ (93 - 3g 2 )piP3 + (21gi q 2 + 6q 3 qj)pi 
+ (7qf + 18gi<?2 + 6939192 - 4gig|)p 3 

They form a quasi bi-inverse-Hamiltonian chain ()6 .3|) 



h (2) - i 
"2 — 3 



9 1 n 2 
UPl ~ fllP3P2 

19391 92 + 6giq2 + T?i)P2 
?g§ - 72g 3 9i92 - 90g?gf - 



(1) 

(2) 
UJ\X\ 

(2) 

Wia; 2 



, _(1) 
P 1 UJqX 1 



P (l>2), , „(2) 



(2) 
U)qX 2 



.(2,1)., _(1) 

1 ^>ox 1 



^(2,2) (2) 

'1 ; w xi 



F 



(2,1) 



(1) 



i 7 



(2,2) (2) 



with the non-canonical symplectic form 



where A = ^p 2 
coefficients 



-,P3 



( 


-B -C 


93 


-A 


-29? \ 


B 


24g? 


-39l 


-3? 2 + 93 





C 


-24g 2 


291 


92 


93 


-<?3 


3qi -2gi 











A 


3<72 - 93 -92 








±9l 


V 2g? 


-93 







/ 


,5 = 


549i92 + 24<j 1 q 3 - 


3pi, C 


= -24gig 2 - 


- 129193 



12qiq 3 + pi and the expansion 



F 



F 



(2,1) _ 



(1,1) 



-93, F, 

„2 1 G„2 



(1,2) 



-91, F 



(2,1) 



,P2 



9?, 



(2,2) 



-2g 3 + 3g 2 , 



- 5939J 1 + 69^92 - 91P1 - 593P2, F 2 (2 ' 2) = -4^ - + 3<? 2 93 + §9iP2 + 9lP3- 
On the extended phase space M. of dimension eight, with additional coordinates ci and C2, the 
differentials dh[ k ^ form a bi-inverse-Hamiltonian chains (|6.7[) 



n Yl 



(i) 



(i) 









niy/ 



ID 



(2) 



(1) 



(2) 
(2) 



with the presymplectic forms 





' w 















v 



d-compatiblc with respect to IIq and vector fields 





dh 
dh 



(2) 


= n lY ^ 


(2) 
2 


= ttiif) 





= 



dh 



ID 



(2) 



(1) 



(o, 0, 1, or, r (2) = (0, 0, 0, iy , = *r - Fr% (i > - pr% 



(i) 



-(I) 



.(1,1)^(1) 



r.(l,2) v (2) 



y(2) = ^(2) 



17,(2, l)v(l) p^/;y 
r l *0 r l 2 C 



(2,2) v (2) 



(2) 



A' 



(2) 



p(2,l)y(l) 



(2,2) v (2) 



-To If 



n 



The bi-Hamiltonian extensions of systems from presented examples the reader can find in j!8) 
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